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CONVERGENCE OF TIME-INHOMOGENEOUS 
GEODESIC RANDOM WALKS AND ITS APPLICATION 
TO COUPLING METHODS 

By Kazumasa Kuwada 1 

Ochanomizu University 

We study an approximation by time-discretized geodesic ran- 
dom walks of a diffusion process associated with a family of time- 
dependent metrics on manifolds. The condition we assume on the 
metrics is a natural time-inhomogeneous extension of lower Ricci cur- 
vature bounds. In particular, it includes the case of backward Ricci 
flow, and no further a priori curvature bound is required. As an ap- 
plication, we construct a coupling by reflection which yields a nice 
estimate of coupling time, and hence a gradient estimate for the as- 
sociated semigroups. 

1. Introduction. It has been well known that there is a strong connection 
between behavior of heat distributions or Brownian motions and geometry 
of their underlying space. Even on time-inhomogeneous spaces such as Ricci 
flow, this guiding principle has been confirmed through recent developments 
(see [1, 8, 17-19, 26, 32] and references therein). Some of them [1, 18] are 
based on coupling methods of stochastic processes. Given two stochastic pro- 
cesses Y\(t) and Y^it) on a state space M, a coupling X(i) = (Xi(t),X2(t)) 
of Y\(t) and is a stochastic process on M x M such that X{ has the 

same law as Yi for i = 1,2. By constructing a suitable coupling which reflects 
the geometry of the underlying structure, one can obtain various estimates 
for heat kernels, harmonic maps, eigenvalues etc. under natural geometric 
assumptions (see [12, 15, 30], e.g.). Since coupling of random variables pro- 
vides a coupling of their distributions, coupling methods are naturally con- 
nected with the theory of optimal transportation, which are used in some 



Received July 2010; revised March 2011. 

1 Supported in part by the JSPS fellowship for research abroad and by KAKENHI 
(22740083). 

AMS 2000 subject classifications. Primary 60F17, 53C21; secondary 58J65, 53C44, 
58J35. 

Key words and phrases. Geodesic random walk, Ricci flow, diffusion process, coupling. 



This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Probability, 
2012, Vol. 40, No. 5, 1945-1979. This reprint differs from the original in 
pagination and typographic detail. 



1 



2 



K. KUWADA 



of aforementioned results [19, 26]. With further studies in this direction in 
mind, here we consider an approximation of diffusion processes associated 
with a family of time-dependent metrics by so-called geodesic random walks. 
Generally speaking, one of the major reasons that we establish approxima- 
tion is to overcome technical difficulties in studying the object in the limit. 
This is also our case, and we will use the approximation in order to study 
a coupling of diffusion processes. 

Let M be a smooth manifold with a family of complete Riemannian 
metrics g(t) indexed by t 6 [Ti,T2]. By (A(i)) tg [ TljT2 ] , we denote the g(t)- 
Brownian motion. It means that X(t) is a time-inhomogeneous diffusion 
process on M associated with A g ^/2, where A g (t) is the Laplacian with re- 
spect to g(t) (see [8] for a construction of ^(t)-Brownian motion). A geodesic 
random walk X a on M with a parameter a is a discrete time Markov 
chain whose one-step variation is given as follows: Given a position x at 
some time t, consider a random vector in T X M. We map it to M by g(t)- 
exponential map to determine the next position. Here the parameter a is 
implemented as a (diffusive) scaling on time step and on the length of the 
random vector in T X M; see Section 3 for more details. In this paper, we con- 
sider only the case that all the random vectors in tangent space is specified 
to the one having a uniform distribution on a <7(t)-ball whose radius is com- 
parable to a. A simplified version of our main theorem, the convergence of 
geodesic random walks, is stated as follows; see Theorem 3.1 and Section 3 
for a more precise and general statement: 

Theorem 1.1. Suppose 
(1.1) d t g(t) < Ric g{t) 

holds. Then a continuous time interpolation of X a converges in law to X as 
a — > 0. 

As we will see in the sequel, there are several technical difficulties arising 
from the time-dependency on the metric. Nevertheless, the assumption of 
the full statement in Theorem 3.1 is much weaker in some respect than that 
in the classical time-homogeneous case. Thus this assertion itself would be 
of interest, independently of its application to coupling methods. 

In the time-homogeneous case, the convergence in law of scaled geodesic 
random walks to the Brownian motion is used to study a coupling of Brown- 
ian motions (Xi(t), A^i)) by reflection; see [16, 27]. A coupling of this kind 
provides us a useful control of the coupling time r* , the first time when X\ 
and X2 meet. Even in our time-inhomogeneous case, Theorem 1.1 carries 
the same estimate in almost the same way. A simplified version of this as- 
sertion is as follows; for the complete statement of our main theorem, see 
Theorem 4.1. 
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Theorem 1.2. Suppose (1.1). Then, for each x±,x 2 G M , there exists 
a coupling X(t) := (X\{t) , X2{t)) of two g(t)-Brownian motions starting at 
(xijXz) satisfying 



(1.2) P[r*>i]< 



inf B(s)> 

Ti<s<t 



for each t, where dg/^A is the distance function on M with respect to g{T\), 
and B(t) is a one- dimensional standard Brownian motion starting at the 
time T\. 

Similarly to the time-homogeneous case, Theorem 1.2 yields a gradient 
estimate of the heat semigroup, which tells us a quantitative estimate on 
the smoothing effect of the heat semigroup; see Corollary 4.3. In addition, 
we can apply our method to construct different kinds of couplings. As one of 
them, coupling by space-time parallel transport is studied in [18] by using 
Theorem 1.1, and it sharpens the result by Topping [26], concerning the 
monotonicity of a transportation cost between the heat distributions whose 
cost is measured by Perelman's J2?-distance. 

Condition (1.1) is essentially the same as backward super Ricci flow in [19]; 
our condition is slightly different in constant since our g(i)-Brownian mo- 
tion, and hence our heat equation corresponds to A g ^/2 instead of A fl ( t ). 
Obviously, (1.1) is satisfied if g(t) evolves according to the backward Ricci 
flow dtg(t) = Ric ff ( t ). From a different point of view, condition (1.1) can be 
interpreted as a time-inhomogeneous analog of nonnegative Ricci curvature 
since dtg(t) vanishes if g(t) is independent of t. Along this viewpoint, we 
can consider a time-inhomogeneous analog of more general lower Ricci cur- 
vature bounds, and we can obtain the conclusion under such generalized 
conditions in the sequel; see Assumption 1 and (4.1); cf. Remark 4.2. It 
should be remarked that, even in those cases, no uniform lower bound of 
Ric 9 ( t ) only in terms of git) without time derivative is assumed. In par- 
ticular, no bounds of g(t) or <7(t)-curvature tensor being uniform in time 
are required. Since a Ricci flow will produce a singularity in a finite time, 
a time-uniform bound on g(t) or Ric ff ( t ) seems to be restrictive. It might be 
possible to simplify the proof by supposing additional assumptions involving 
a time-uniform estimate; however, this is out of the scope of this paper. 

In our argument, the distance function d g u\ with respect to the time- 
dependent metric g(t) plays a prominent role. The first variation in t of d g n\ 
is described in terms of dtg(t), and the second variation of d g u\ in space 
variables involves a notion of curvature. Both of these variations appear in 
the bounded variation part of the radial process d g u\{o, X(t)) of the g(t)- 
Brownian motion via the ltd formula. Thus a relation between dtg(t) and 
Ric s ( t ) [e.g., (1.1)], produces a nice control of the radial process. Although 
we will work on geodesic random walks instead of the g(i)-Brownian motion 
itself, such an observation is still efficient. 
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In the time-homogeneous case, the convergence of scaled geodesic ran- 
dom walks is first studied by J0rgensen [14] by using the convergence theory 
of semigroups; see [6, 22], also. However, in our framework, it is not clear 
whether we can apply a similar technique since the base measure, the Rie- 
mannian volume, depends on time, and hence we cannot expect that it makes 
the heat semigroup invariant. To avoid such a technical difficulty, we use the 
uniqueness of the martingale problem instead for identifying the limit. An- 
other difficulty arises from the lack of time-uniform bounds of Riemannian 
metrics. It prevents us to expect a global comparison of geometric structures, 
such as d g M, between different times. Thus we will make some efforts for 
localizing the problem by giving a uniform estimate of the first exit time 
of X a from a large ball centered at a reference point. Note that our assump- 
tion admits lower unbounded Ricci curvatures even in the case dtgit) = 
(see Assumption 1). Thus our assumption on the geometry of the underlying 
space is weaker than that in [14] (by considering Riemannian manifolds with 
lower unbounded curvature, we can easily find an example which does not 
satisfy the assumption in [14]). On the other hand, the assumption on the 
driving noises of the geodesic random walk in [14] is more general than our 
specified one. Though it might be possible to take a more general noise under 
our assumption, our result already works well for applying the approxima- 
tion to coupling methods. As a related work, the theory of time-dependent 
Dirichlet forms has been developed for studying the time-inhomogeneous 
Markov processes in the literature [20]; see also [24]. Unfortunately, because 
of above-mentioned difficulties, our framework does not fall into the scope of 
those theories at this moment. It might be an important problem to extend 
those theories so that they includes our case. 

The organization of this paper is as follows. In the rest of this section 
after this paragraph, we review existing approaches on the construction of 
couplings. By comparing those approaches with ours, we try to explain the 
reason why we choose our approach for constructing a coupling by reflec- 
tion. In the next section, we show basic properties of a family of Rieman- 
nian manifolds ((M, g{t)))t- In particular, we prove that Riemannian metrics 
(g(t))t are locally comparable with each other. It will be used to give a uni- 
form control of several error terms which appear as a result of our discrete 
approximation. In Section 3, we will study geodesic random walks in our 
time-inhomogeneous framework. There we introduce them and prove the 
convergence in law to a diffusion process. After a small discussion at the be- 
ginning of the section, the proof is divided into two main parts. In the first 
part, we will give a uniform estimate for the exit time of geodesic random 
walks from a big compact set. Our assumption here is almost the same as 
in [17] where nonexplosion of the diffusion process is studied; see Remark 3.3 
for more details. In the second part, we prove the tightness of geodesic ran- 
dom walks on the basis of the result in the first part. In Section 4, we will 
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construct a coupling by reflection and show an estimate of coupling time, 
which completes the proof of Theorem 1.2 as a special case. In Section 5, we 
will give a short remark about how our method is also applicable to study 
a coupling by parallel transport. 

1.1. Existing arguments on coupling methods. As stated above, we com- 
pare our method of the proof with existing arguments in coupling methods 
from a technical point of view. We hope that the following observation will 
be helpful to extend coupling arguments other than our own in this time- 
inhomogeneous case. 

In order to go into details, let us review a heuristic (and common) idea of 
the construction of a coupling by reflection as well as that of the derivation 
of (1.2). Given a Brownian particle X±, we will construct X2 by determining 
its infinitesimal motion dX2(t) G Ty 2 (t)M by using dX\(t) G T^wM. First 
we take a minimal g(i)-geodesic 7 joining X\(t) and X2{t). Next, by using 
the parallel transport along 7 associated with the g(t)-Levi-Civita connec- 
tion, we bring dX\{t) into Tx 2 (t)M. Finally we define dX-i{t) as a reflection 
of it with respect to a hyperplane being g(i)-perpendicular to 7 in T X2 (t)M. 
From this construction, the ltd formula implies that d g ^(Xi{t),X2(t)) should 
become a semimartingale at least until (Xi(t),X2(t)) hits the #(i)-cutlocus 
CutgM. The semimartingale decomposition is given by variational formu- 
las of arc length. On the bounded variation part, there appear the time- 
derivative of d g u\ and (a trace of) the second variation of d g u\, which is dom- 
inated in terms of the Ricci curvature. With the aid of our condition (1.1), 
these two terms are compensated and a nice domination of the bounded vari- 
ation part follows. Thus the hitting time to of d g ^(Xi(t),X2(t)), which 
is the same as r* , can be estimated by that of the dominating semimartin- 
gale. Indeed, we can regard 2B{t) + d g (T 1 \(xi,X2) which appeared in the 
right-hand side of (1.2) as the dominating semimartingale. The effect of 
our reflection appears in the martingale part 2B(t) which makes it possi- 
ble for the dominating semimartingale to hit 0. This construction seems 
to work as long as (Xi(t), A^(i)) is not in the cutlocus. Moreover, if we 
succeed in constructing it beyond the cutlocus, then the same domination 
should hold. Indeed, the effect of singularity at the cutlocus should decrease 
d g ^(Xi(t),X2(t)). Thus a "local time at the cutlocus" will be nonpositive, 
and hence negligible. 

After this observation, we can conclude that almost all technical diffi- 
culties are concentrated on the treatment of singularity at the cutlocus in 
order to make the heuristic argument rigorous. In fact, Theorem 1.2 is shown 
in [21] by using SDE methods under the assumption that the g(t)-cutlocus is 
empty for every t € [Ti,T2]. It should be remarked that the joint distribution 
of the coupled particle (Xi(i), A^t)) could be singular to the Riemannian 
measure onMxM (at least it is the case when M is a flat Euclidean space) . 
Thus it is not clear that the cutlocus is really "small" for the coupled par- 
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tide despite the fact that the cutlocus (as a subset of M x M) has null 
g(i)-Riemannian measure. 

In our approach, we first construct a coupling of geodesic random walks 
and then take a limit to obtain the desired coupling. Since we first derive 
a dominating semimartingale for coupled geodesic random walks, we need 
only a difference inequality instead of the Ito formula. By virtue of this 
difference, we can obtain a desired estimate beyond the cutlocus by dividing 
a minimal geodesic joining particles into small pieces so that the endpoints 
of each piece are uniformly away from the cutlocus; see Lemma 4.4. As 
a result, we can avoid extracting a local time at the cutlocus and directly 
obtain a dominating process which does not involve such a term. Moreover, 
the dependency on time parameter of the cutlocus does not cause much 
difficulty in our approach. 

In the time-homogeneous case, there are several arguments [9, 12, 28-30] 
to construct a coupling by reflection by approximating it with ones which 
move as mentioned above, if they are distant from the cutlocus and move 
independently if they are close to the cutlocus. In some of those arguments, 
we need to estimate the size of the total time when particles are close to 
the cutlocus. In such a case, an extension of these arguments to the time- 
inhomogeneous case does not seem straightforward since the </(£)-cutlocus 
depends on time and estimates should be more complicated. The argument 
in [30] uses supermartingales to extract the local time at the cutlocus in an 
implicit way, and no estimate of times spent around the cutlocus is necessary. 
Thus it seems possible to extend his argument in the time-inhomogeneous 
case. Since his argument relies on some detailed properties of parabolic 
PDEs, we need to develop time-inhomogeneous analogs of them to com- 
plete this plan. The fact that our assumption (1.1) [or (4.1)] does not imply 
any time-uniform lower bound of the Ricci curvature by a constant might 
be an obstacle. 

If we employ the theory of optimal transportation, we will work on cou- 
plings of heat distributions instead of coupling of Brownian motions. Once 
we move to the world of heat distributions, we can expect that the cutlo- 
cus is treated more easily since they are of measure zero with respect to 
the Riemannian measure. However, at this moment, the theory of optimal 
transport is not so strong a tool in this context for the following two rea- 
sons. First, the range of the theory is restrictive in the sense that it only 
deals with couplings corresponding to the coupling by parallel transport. 
Second, the theory of optimal transportation provides a weaker result than 
a probabilistic approach does, even in studying couplings by parallel trans- 
port; for instance, see [19] and compare it with [1]. It should be remarked 
that such a difference between these two approaches exists even in the time- 
homogeneous case. 

Arnaudon, Coulibaly and Thalmaier [1] recently developed a new method 
to construct a coupling, which works even in the time-inhomogeneous case. 
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They consider a one-parameter family of coupled particles along a curve. 
Intuitively speaking, they concatenate coupled particles along a curve by 
iteration of making a coupling by parallel transport. Since "adjacent" par- 
ticles are infinitesimally close to each other, we can ignore singularities on 
the cutlocus when we construct a coupled particle from an "adjacent" one. 
It should be noted that their method does not seem to be able to be ap- 
plied directly in order to construct a coupling by reflection. Indeed, their 
construction of a chain of coupled particles heavily relies on a multiplicative 
(or semigroup) property of the parallel transport. However, our reflection 
operation obviously fails to possess such a multiplicative property. Since our 
reflection map changes orientation, there is no chance to interpolate it with 
a continuous family of isometries. 

2. Properties on time-dependent metric. As in Section 1, let M be a Tri- 
dimensional manifold and (g(t))te[T 1 ,T 2 ] a family of complete Riemannian 
metrics on M which smoothly depends on t, for — oo < T\ < T 2 < 00. 

Remark 2.1. It seems to be restrictive that our time parameter only 
runs over the compact interval [7\,T2]. An example of g(t) we have in mind 
is a solution to the backward Ricci flow equation. In this case, we can work 
on a semi-infinite interval [Ti,oo) only when we study an ancient solution 
of the Ricci flow. Thus T2 < 00 is not so restrictive. In addition, we could 
extend our results to the case on [Ti,oo) with a small modification of our 
arguments. It would be helpful to study an ancient solution. To deal with 
a singularity of Ricci flow, it could be nice to work on a semi-open interval 
(Ti,?^]) where T\ is the first time when a singularity emerges. In that case, 
we should be more careful since we cannot give "an initial condition at Ti" 
to define a (?(/;) -Brownian motion on M. 

We collect some notation which will be used in the sequel. Throughout 
this paper, we fix a reference point o G M. Let No be nonnegative integers. 
For a,i)Gffi, a A 6 and a V b stand for min{a, 6} and max{a, &}, respectively. 
Let Cut s ( t )(x) be the set of the <?(i)-cutlocus of x on M. Similarly, the 
<?(£) -cutlocus Cut 9 (j) and the space-time cutlocus Cutsx are defined by 

Cut 9(t) := {(x,y) G M x M\y G Cut g{t) (x)}, 

CutgT :={(t,x,y) G [T U T 2 ] xMx M\(x,y) G Cut g(i) }. 

Set D(M) := {(x,x)\x G M}. The distance function with respect to g(t) is 
denoted by d g ^(x,y). Note that CutsT is closed and that d g ^(-, •) is smooth 
on [Ti, T 2 ] x M xM\ (Cut ST U [T u T 2 ] x D(M )); see [19]; cf. [17]. We denote 
an open ^(s)-ball of radius R centered at x G M by Bj^\x). Some additional 
notation will be given at the beginning of the next section. 
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In the following three lemmas (Lemmas 2.2—2.4), we discuss a local com- 
parison between d g u\ and d g ^ for s 7^ t. Those will be a geometric basis of 
further arguments. 

Lemma 2.2. Let Mq be a compact subset of M. Then there exists k = 
k(Mq) such that 

e- 2K ^g(s)<g(t)<e 2K ^g( S ) 

holds on Mq for t,s G [Ti,T2]. In particular, if a minimal g(s)-geodesic 7 
joining x,y G Mq is included in Mq, then, for t G [T\,T2\, 

d g{t) (x,y) < e^-^dg^x^). 
Proof. Let ir:TM — > M be a canonical projection. Let us define Mq 

by 

M :={(t,v) G [Ti,T 2 ] x TM\tt{v) G M ,\v\g(t) < !}■ 

Note that Mo is closed since g(-) is continuous. We claim that Mq is sequen- 
tially compact. Let us take a sequence ((t n ,u n )) nS M C Mq. We may assume 
£ n — >■ i G [Ti,T2] and 7r(u n ) Mo as n -> 00 by taking a subsequence if 

necessary. Let [/ be a neighborhood of p such that {v G TM|7r(-u) G £/} ~ 
U x M m . For sufficiently large n, v n is in U x M m and we write v n = (p n ,v n ). 
If we cannot take any convergent subsequence of (f n )neN> then \v n \ — > 00 as 
n — > 00, where | • | stands for the standard Euclidean norm on R m [irrelevant 
to (s(i))te[Ti,T 2 ]]- Set vj, = (p n , |v n | _1 ^n)- Then, there exists a subsequence 
( u n fc )fceN C (v^)ngN such that u' — >■ u^j = (p, «') as n — > 00 for some u' G M m 
with \v'\ = 1. Since <?(•) is continuous, g(t nk )(v' ,v' ) — > ffCOOfc,^) as — >• 
00. On the other hand, g(t nk )(v' nk ,v' nk ) < \v nk \~ 2 -> since g(t n )(v n ,v n ) < 1. 

Thus u' must be 0. It contradicts with = 1. Hence Mo is sequentially 
compact. 

Since Mo 9 (i, v) h-> dtg(t)(v,v) is continuous, there exists a constant k = 
k(M ) >0suchthat |%(t)(u,v)| < 2k for every (t,u) G M . Take u G 7r -1 (M ), 
v^0 w{v) . Then 

d t g(t)(v,v) = \v\ 2 g{t) dtg{t){\v\-l- t) v,\v\-^ t) v) < 2 K \v\ 2 g{t) . 

Thus dt\ogg(t)(v,v) < 2k holds. By integrating it from s to i with s <t, 
we obtain g(t)(v,v) < e 2 *^^ g(s)(v,v). We can obtain the other inequality 
similarly. 

For the latter assertion, for a, 6 with 7(a) = x and 7(6) = y, 

rb rb 

d g{t) (x,y)< |7(u)| 5(t) dn<e K l*- s l / \j(u)\ g(s) du 



■ e K| *" 
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Lemma 2.3. For R>0, x G M and t G [T 1 ,T 2 ], there exists 5 = 6(x,t, 
R)>0 such that } (x) C B$(x) forr<R and s G [Ti,T 2 ] with \s-t\< 5. 

PROOF. Set k := «(B^(a;)) as in Lemma 2.2 and 5 := « -1 log2. Take 
p G Br(x) and a minimal g(s)-geodesic 7: [a, b] — > M joining x and p. Sup- 
pose that there exists no G [a, 6] such that 7(^0) G B^! (x) c . Let uo := inf{u € 

[a,6]| 7 (u) G ^(rc) 6 }. Since 7 ([«,«o]) C B®(x) C B§(x) and d g(t) (x, 
7(7x0)) = 3r, Lemma 2.2 yields 

'"• • - - KS r° . . 3r 



dg( s )(x,p)> \j(u)\ g ^du>e K / \i{u)\ g (t)du> — 

J a J a 

This is absurd. Hence 7([a,6]) G B^(x). In particular, 7(6) = p£ B^(x). 
□ 

Lemma 2.4. For i? > 0, f/iere exists a compact subset Mq = Mq(R) of M 
such that 

(2.1) \ P eM\ inf d fl(t) (o,p)<J?}cM . 

Proof. For each i G [Ti,T 2 ], take 5(o,t,i?+ 1) > according to Lem- 
ma 2.3. Take {U}f =1 C [T U T 2 ] such that 

n 

[T U T 2 ] C Qfc- 5(o,ti,R + l), U + 6(o,U,R + l)). 
i=i 

Let us define a compact set Af C M by M := |J"=i B { ^{o). Take p G M 
such that infr 1 <i<r 2 d g ^(o,p) < R. For e G (0, 1), take s G [Ti,T2] such that 
dg( s ){°iP) < R + £■ Then there exists j G {1, . .. ,N} such that \s — tj\ < 

5(o,tj,R + 1). By Lemma 2.3, it implies p G B^]_ e (o) C bM +£ Jo) C 
UiLi -^3(B+ e )(°)- Hence the conclusion follows by letting £4.0. □ 

Another useful consequence of Lemmas 2.2 and 2.3 is the following: 

Lemma 2.5. d g r.\(-,-) is continuous on [Ti,T 2 ] x M x M. 

Proof. Since the topology on [7\,T2] x M x M is metrizable, it suf- 
fices to show lim^oo dg^) (x„, y n ) = d g ^)(x,y) when (t ) -> (i,x,y) 
as n — )■ 00. By the triangle inequality, 

K(t n )(>n,yn) -d<?( t )(x,y)| < |d ff ( tn )(x,y) - 

(2.2) 

+ d g{tn) (x,x n ) + d g(tn) (y,y n ). 
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Take R > so that (x) includes a minimal (j(i)-geodesic joining x and y. 

Take k = k(B^(x)) according to Lemma 2.2. We can easily see that every 

minimal g(i)-geodesic joining y and y n is included in B^ix) for sufficiently 
large n G N. Thus Lemma 2.2 yields 

limsupd ff ( tn )(y,y n ) < limsupe K| *" in| (i 9 ( t ) (y, y n ) = 0. 

n— >oo n— >oo 

We can show d g u\(x, x n ) —> similarly. Take a minimal g r (t n )-geodesic 
7 n : [a, 6] — > M joining x and y. By our choice of R, Lemma 2.2 again yields 

d g (t n )(x,7n(u)) < d g(tn) (x,y) < e K|t -*" l d 9(t) (x,y) < e^'^R. 

It implies limsup n _ i , 00 (i 9 ( 4n )(2;,y) < d g ^(x,y). In addition, j n is included in 

B^j 3 (x) for sufficiently large n. Thus Lemmas 2.3 and 2.2 yield d g ^(x,y) < 

e K '* _ *"'d g ( tn )(x,y). Hence the conclusion follows by combining these esti- 
mates with (2.2). □ 

Before closing this section, we will provide a local lower bound of the 
injectivity radius which is uniform in time parameter. 

Lemma 2.6. For every M\<ZM compact, there is fo=fo(Mi) > such 
that d g (£) (y,z)< f o implies (t, y, z) £ CutsT for any (t, y, z) G [Ti , T%] x M\ x M\ . 

PROOF. Take R > 1 so that sup tg [ Tl r i sup^g^ d g ^(o,x) < R — 1. By 
Lemma 2.4, there exists a compact set Mo C M such that (2.1) holds. For 
every t G [T\,T2] and x G Mi, (t,x,x) ^ CutsT- It implies that there is rj t . x G 
(0, 1) such that (s,y,z) £ CutsT whenever d g ^(x,y) V d g ^(x,z) V |t — s| < 
r/ 4jX since Cutgx is closed. Thus there exist N G N and (ti, Xi) G [Ti, Tq\ x Mi 
(i = 1, . . . , iV) such that 



N , v 

[Ti,T 3 ] x Mi C U - ^p, ^ + x B 



Set f o > by 



1 / k 



rn := — exp max r». T . mm 77/. r ., 

u 2 P V 2i<i<Jv' t " I, /i<Kiv"" 1 " 



where k = k(M ) > is as in Lemma 2.2. Take (s, y, z) G [Ti, T 2 ] x Mi x Mi 
with d g{s )(y,z) < f . Take j G {l,...,iV} so that |s - V d g ( tj) (xj,y) < 
Vtj,xj/2- By virtue of the choice of R and Mo, Lemma 2.4 yields that every 
g(s)-geodesic joining y and 2 is included in Mo. Thus Lemma 2.2 yields 

d gitj) (y,z) <e^-^d g(s) (y,z) < 

It implies \s - tj\ V d g{tj )(xj,y) V d g (t.)(xj,z) < r] t])Xj and hence (s,y,z) ^ 

CutsT- □ 
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3. Approximation via geodesic random walks. Let (-Z(i))te[Ti,T 2 ] be a fam- 
ily of smooth vector fields continuously depending on the parameter t G 
[Ti , T 2 ] . Let X(t) be the diffusion process associated with the time-dependent 
generator Jz? t = A g r t \/2 + Z(t); see [8] for a construction of X(t) by solving 
a SDE on the frame bundle. Note that (t,X(t)) is a unique solution to 
the martingale problem associated with dt + Jz^ on [Ti,T2] x M; see [12] 
for the time-homogeneous case. Its extension to time-inhomogeneous case is 
straightforward; see [25] also. 

In what follows, we will use several notions in Riemannian geometry such 
as exponential map exp, Levi-Civita connection V, Ricci curvature Ric etc. 
To clarify the dependency on the metric g(t), we put (t) on superscript 
or g(t) on subscript. For instance, we use the following symbols: exp"', 
and Ric g (£) . We refer to [7] for basics in Riemannian geometry which will be 
used in this paper. 

For each t G [Ti,T 2 ], we fix a measurable section :M ->• G^{M) of 
the <7(t)-orthonormal frame bundle G^(M) of M. Take a sequence of in- 
dependent, identically distributed random variables {£ n }ngN which are uni- 
formly distributed on the unit disk in W n . Given xq G M, let us define 
a continuously-interpolated geodesic random walk (-X" a (i))te[Ti,T 2 ] on M 
starting from xq with a scale parameter a > inductively. Let tn := (Ti + 
a 2 n) A T 2 for n G N . For t = T 1= , set X a (Ti) := x . After A Q (t) is 
defined for t G [Ti,t^ ], we extend it to t G [t^ j^i+i] by 

| n+1 := V^T2^ (e)) (^(4 a) ))en+i, 

X«(t) := exp^ )} (iZ^ofc*.! + a 2 Z(4«)))) . 

For later use, we define iV (a) := inf{n G N |4+i ~ ^ < a 2 }. This is the 
total number of discrete steps of our geodesic random walks with scale pa- 
rameter a. Set <T := C{[T 1 ,T 2 ] -»• M) and := D([T 1 ,T 2 ] -> Af), the space 
of right continuous paths on M parametrized with [7\ , T 2 ] possessing a left 
limit at every point. By using a distance dgm) on M, we metrize & and ^ 
as usual so that ^ and ^ become Polish spaces; see [10] for a distance 
function on @, for example. Set := C([T\,T 2 ] — > [0,oo)). Let us define 
a time-dependent (0, 2)-tensor field (VZ(t)y by 

(VZ(t)f(X,Y) := j((vWz(t),y) fl(t) + (vg>Z(t), A) 9(i) ). 

Assumption 1. There exists a locally bounded nonnegative measurable 
function b on [0, oo) such that: 

(i) For allfG [T X ,T 2 ), 

2(VZ(i)) b + d t g(t) < Ric g{t) +b(d g(t) (o, -))g(t). 
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(ii) For each C > 0, a one-dimensional diffusion process yt given by 



dy t =dp t + Uc+ jT b(s) ds\ dt, 



where fit is a standard Brownian motion, does not explode. (This is the case 
if and only if 

J exp^— J h{z)dz\ j exp^y b(£) d$\ dzdy = oo, 
where b(y) := C + f^b(s)ds; see, e.g., [13], Theorem VI.3.2.) 

Note that (1.1) is a special case of Assumption 1. Now, we are in position 
to state the main theorem of this paper. 

Theorem 3.1. Under Assumption 1, X a converges in law to X in 
as a — > 0. 

Most of arguments in this section will be devoted to show the tightness, 
that is: 

Proposition 3.2. (^ a ) Qe (o,i) is tight in . 

In fact, as we will see in the following, Proposition 3.2 easily implies 
Theorem 3.1. 

Proof of Theorem 3.1. By virtue of Proposition 3.2, for any sub- 
sequence of (^ a ) a e(o,i) there exists a further subsequence (X° fe ) fcgN which 
converges in law in as k — > oo. Thus it suffices to show that this limit has 
the same law as X. Let (/3 a (t))te[o,oo) De a Poisson process of intensity a~ 2 
which is independent of {£ n }neN- Set 

:={T x + o?^{t-T x ))M%. 

Then the Poisson subordination X ak (/3 ak (-)) also converges in law in @ to 
the same limit; see [5], for instance. Note that (/3 Q (i), ^ Q (/3 a (£)))te[Ti,T 2 ] is 

a time-inhomogeneous Markov process. The associated semigroup P^ and 
its generator Jz? ^ are given by 



P («) 



Z>iV( Q ) 



J?^f:=a- 2 (q^f-f), 
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where 

q {a) f(t, x) := E[f(t + o^,exp^\ay/m + 2^(x)^ 1 + a 2 Z(t)))}. 

We can easily prove Jzf ^ a 'f — > (dt + »5f )/ uniformly as a — > for / G Cg°([Ti, 
T 2 ] x M). Since (/3 a (i), X Q (/3 a (£))) tg[TliT2] is a solution to the martingale 
problem associated with J£^ a \ the limit in law of (fi ak {■) , X a (J3 ak (■))) solves 
the martingale problem associated with dt + J£.. By the uniqueness of the 
martingale problem, this limit has the same law as that of (i, ^(i))te[Ti,T 2 ] • 
It completes the proof. □ 

Remark 3.3. Proposition 3.2 also asserts that any subsequential limit 
in law is a probability measure on c f. Since we have not added any cemetery 
point to M in the definition of Theorem 3.1 implies that X cannot 
explode. It almost recovers the result in [17]. Our assumption is slightly 
stronger than that in [17] on the point where we require (ii) for all C > 0, not 
a given constant. Note that we will use Assumption 1 (ii) only for a specified 
constant 2Cq given in Lemma 3.9. However, its expression looks complicated, 
and it seems to be less interesting to provide an explicit bound. 

Now we introduce some additional notation which will be used in the rest 
of this paper. For t G [Ti,T 2 ], we define [t\ a by 

[t\ a := sup{a 2 n + T\\n G No, a 2 n + T\ < t}. 

Set & n := cr(£i, . . . ,£ n ). For R>1, let us define a R : [T Xl T 2 ] U {oo} by 

a R (w) := inf{t €[T 1 ,T 2 ]\w(t)>R-l}, 

where inf = oo. We write a R := a R (d g ^(o, X a (-))) and <tr := a~ 2 {[a R \ a — 
T\) + 1. Note that a R is an J?" n -stopping time. For each t G [Ti,T2] and x, 

y G M with x / y, we choose a minimal unit-speed g(i)-geodesic j X y '■ [0, 

d g (t)(x,y)] M from x to y. Note that we can choose jxy so that (x,y) \— > 

jxy is measurable in an appropriate sense; see, for example, [27] . We use the 

same symbol 7^ for its range ^/xy{[0,d g ^(x,y)]). 

3.1. A uniform bound for the escape probability. The goal of this subsec- 
tion is to show the following: 

Proposition 3.4. lim Rtoo limsup Q ^ P[<7R < T 2 ] =0. 

For the proof, we will establish a discrete analog of a comparison argu- 
ment for the radial process as discussed in [17]. From now on, we fix R > 1 
sufficiently large so that g? 9 (Ti) (°i ^0) < -R — 1 until the final line of the proof 
of Proposition 3.4. We also fix a compact set Mq C M satisfying (2.1). Set 
r o := ^0 A (1/2), where fo = fo(Mo) is as in Lemma 2.6. 
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The first step for proving Proposition 3.4 is to show a difference inequality 
for the radial process d g r t \(o, X a (t)) (Lemma 3.7). It will play the role of 
the ltd formula for the radial process in our discrete setting. We introduce 
some notation to discuss how to avoid the singularity of d g ^(o, •) on {o} U 
Cut s (.)(o). For r > 0, let us define a set A' r ,A" and A r as follows: 

A' r :={(t,x,y)e[T 1 ,T 2 ]xM x M \d g{t) (x,x') + d g(t) (y, y') + \t - 1'\ > r 

for any (t',x',y') G Cut ST }, 

A';:={(t,x,y)e[T 1 ,T 2 ]xMoxM \d g{t) (x,y)>r}, 

Note that A r is compact and that d g ^(-, •) is smooth on A r . For t G [Ti,T2] 
and p G M, let us define Op ^ G Mq by 

(t) := J 7$ (y) . if (*,o,p)^^, 
^ o, otherwise. 

(t (0!) ) 

For simplicity of notation, we denote o n {a) by o n . Similarly, we use the 

(t ) (t) 
symbol j n for 7 {a) throughout this section. Note that (t,o p ,p) f. 

CutsT holds. Furthermore, it is uniformly separated from CutsT in the fol- 
lowing sense: 

Lemma 3.5. There exist r% > and 5\ > such that the following holds: 

let to,te [Ti,T 2 ] with t - t G [0,<*i]. Lei p £ B&i(o) and p G ^ o) (po)- 
T/ien we Ziaue: 

(i) d g(t) (o,p) < e< t ~ tQ \d g{to) (o,p ) + d g{tQ) (po,p)); 

(ii) (t,Op °\p) G A n rtenpo^^^o). 

//ere k = k(Mo) > is given according to Lemma 2.2. 

By applying Lemma 3.5 to X a , we obtain the following: 

Corollary 3.6. There exist ao > and h : [0, ao] — > [0, 1] iwzf/i lim Q ^o h(a) = 

suc/i i/iaf the following holds: for a < ao, n G No and s,t G [4f> w/ien 
n < o^: 

(i) d 9{t) ( ,X a ( S ))<e^ 2 (d 5(i ( Q)) ( ,X«(ti a) )) + / l (a)); 

(ii) (t,o n ,X a (s)) G ^ w/ien X a (t ( n a) ) <£ B^p\o). 
Here r\ is the same as in Lemma 3.5. 
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Proof. Set Z '■= sup t& ^ T ^ xeMo \Z(t)\ g ^(x). Note that we have 

d ,A a) .(X a (t^),X a (t))<V^+2a + Za 2 

by the definition of X a . Take oq > so that \Jm + 2oq + Za^ < 5± and 
a 2 < 5± hold, where <5i is as in Lemma 3.5. Then the conclusion follows by 
applying Lemma 3.5 with to = tn \ Po = X a (tn) and p = X a (s). □ 

Proof of Lemma 3.5. We show that (i) holds with <5i = 1. By the 
triangle inequality, the proof is reduced to showing the following two in- 
equalities: 

(3-1) d g{t) (o, Po )<e^ t - to U gito) (o, Po ); 
(3.2) dg(t)(po,P) < e K{t - to) d g{to) (po,p). 

Our condition (2.1) yields that 7op^ ^ s included in Mq. Thus Lemma 2.2 
yields (3.1). When p G B^°\po), we have yjO$ C B^°\o). Hence (2.1) and 
Lemma 2.2 yield (3.2) in a similar way as (3.1). 



Let us consider (ii). For simplicity of notation, we denote Op^ by d in 
this proof. We assume that t — to G [0,5] and p G Bf°\po) hold for S > 0. 
First we will show (t,o',p) G A" Q , 4 wnen ^ * s sufficiently small. Note that 

(t ,d,p ) G A'^ /2 holds since p £ B?°\o) and d fl ( to )(o,o') G {r /2,0}. Let 

5 ^ To'po ' ^ ^ e triangle inequality, 

(3-3) d g(t) (o,q) < d g{ t){o,d) + d g(t) (d ,p Q ). 

Since r /2 < 1 < ii holds, (2.1) yields C M when d / o. We can easily 
see that 7^°^ C 7^,° C Mq. Thus, by applying Lemma 2.2 to (3.3), 

/ d g(t) (o,q)<e^ t - to \d g{to) (o,d) + d g{to) (d,p )) 

3-4 . 

< (R-l)e KS . 

Take 5 2 := 1 A (k" 1 \og(R/(R — 1))). Then, for any 5 G (0,<5 2 ), (3.4) and (2.1) 
imply 7q/L C Mq. Hence the triangle inequality, Lemma 2.2 and (3.2) yield 

d 9 (t){o',p) > d g( t)(d,po) - d g{t) (p ,p) 

(3.5) > e- K ^d g{to) (d,p ) - e^-^d g{to) (p ,p) 

when 5 < 62- Thus there exists £3 = 6s(K,ro,R) G (0, #2] such that the right- 
hand side of (3.5) is greater than ro/4 whenever 5 G (0,^3). Hence (t,o',p) G 
A" ,. holds in such a case. 

ro/4 



16 



K. KUWADA 



Next we will show that there exists r[ > such that (t,o',p) G A ' , holds 
for sufficiently small 5. Once we have shown it, the conclusion holds with r% = 
r[ A (ro/4). As we did in showing (t, o' ,p) G A" , 4 , we begin with studying the 

corresponding statement for (to,o',po). More precisely, we claim that there 
exists r'l G (0,1) such that (to,o',po) G A r ». When d = o, (to,o',po) G A' ro 

directly follows from the definition of o' = Opf 1 . When o' ^ o, set 

H := {(t,x,y) G [Ti,T 2 ] x M x M |r < d fl(t) (o,j/) <R-1, 

d g(t)(o,x) = r /2,d g ( t )(x,y) = d g{t) {o,y) - d g{t) (o,x)}. 

Note that H is compact and that H n CutsT = holds since (t,x,y) G ii 
implies that x is on a minimal g(i)-geodesic from y to o. Since (to, o',po) G -ff 
by the definition of o', it suffices to show that there exists f\ > such 
that H C .A- . Indeed, the claim will be shown with r'{ = f\ A ro once we 
have proved it. Suppose that H C A' r does not hold for any r G (0, 1). Then 
there are sequences (tj,Xj,yj) G H, (t'j,x'j,y'j) G CutsT, 3 £ N, such that 
|ij - iy| + dg^.jixj,^) + dg^^yj^y'j) ->• as j -> oo. We may assume that 
((tj,Xj,yj))j converges. Since (tj,Xj,yj) G -ff, £j,2/j £ Mo holds for suffi- 
ciently large j. Thus we can take a convergent subsequence of ((t'px'-,y'-))j. 
Since CutsT and if are closed, and d g ^(-,-) is continuous, it contradicts 
with Hf] Cut S T = 0- 

To complete the proof, we show that there exists Si G (0, S3] such that 
(t,o',p) G ^4' r /// 2 when <5 G (0, #i). Suppose that there exists (t',x',y') G CutsT 

such that \t - t'\ + d g (t)(o',x') + d g ^(p,y') < r"/2. For any g G 7®,, the 
triangle inequality and assertion (i) yield 

(3.6) d g(t) (o, q) < d g(t) (o, p) + d m (p, y') <e K5 (R-l + S)+ r'{/2. 

A similar observation implies d g u)(o,q') < (e K< Vo +r")/2 for q' G %) x f Thus 
there is S4 = S^(k, R) G (0,53] such that the right-hand side of (3.6) is less 

than R and (e K< Vo + r'{)/2 < R whenever 5 G (0, 84). In such a case, 7^* , C M$ 
and j% C M hold. Since (to,o',po) gA^,„, Lemma 2.2 yields 

\t-t'\+ dgQ (of, x') + d g (t) (p, y) 

(3.7) > \t -t'\-5 + e- KS d g{to) (o',x') + e- KS d g{to) (p,y') 

> e~ K8 r'{ + (1 - e~ KS )\t -t'\-8- e~ K5 S. 

Take <5i = Si(k, r'{) G (0,^4] so that the right-hand side of (3.7) is greater 
than r'l/2 when 5 G (0,<5i). Then (3.7) is absurd for any 5 G (0,<5i). Thus it 
implies the conclusion. □ 

We prepare some notation for the second variation formula for the arc 
length. Let be the g(i)-Levi-Civita connection and the g(t)-curva- 
ture tensor associated with V^. For a smooth curve 7 and smooth vector 
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fields U,V along 7, the index form Ij\u,V) is given by 

i®(u,v) ■.= [ ({vVu,v®v) g{t) - {nV{u,^,v) g{t) )ds. 

We write Ij (U, U) =: Ij(U) for simplicity of notation. Let Gt :X ,y(u) be the 
solution to the following initial value problem on [0,d(x,y)]: 

G t ,x,y{ u ) ~ nT^l G t,x,y{U), 

<Wo) = o, <r;;.,, v (0) 1. 

Note that Gt, x ,y{u) > for u G (0, d(x, y)] if y ^ Cut 5 ( t )(a;); see [17], proof of 
Lemma 9. For simplicity, we write G n := G j a ) va ,A a h- When X a (i„ ) ^ 
Sr n (o), we define a vector field V* along 7 n for each V G T < a ),M by 



:= — ^ V)( 7n (n)), 



(t^) 

where //\" V is the parallel vector field along 7 n of V associated with 

V ( *"" ). Take v G M. m . By using these notations, for n G No with n < , 
let us define A n+ i and A n+ i by 

A n+ i := d t d ^^0,0^ + d t d g{t ^ )) {o n ,X a {t ( £ ) )) 

+ (Z(4«)),7„) 9(i ( Q)) (^(t(r))) + ) )(4 fl ). 



when X Q (i}T ) ^ -Br " 0), and A n+ i := x/m+^^n+i , u) R m and A n+1 := 
otherwise. 

Lemma 3.7. // n < or A N^ a \ a < «o is sma£Z enough and X a (tn) ^ 
5} n j (o), then 

d a(t W ) (o,^ a (tS 1 )) < ^ ( ^))(o^ a (4 Q) )) + «A n+1 + a 2 A n+1 + o(a 2 ) 

almost surely, where Qo is as in Corollary 3.6. In addition, o(a 2 ) is controlled 
uniformly. 

Proof. By virtue of Corollary 3.6, for sufficiently small a, the Taylor 
expansion together with the second variation formula yields 
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(3.8) 



d 9(t ( Q)) ( 0n ,A^^ 



+ o(a 2 ), 



where i is a ^(tn^-Jacobi field along j n with a boundary value con- 
dition +1 (o n ) = and ^(X a (t^)) = £ n +i- Note that o(a 2 ) can be 
chosen uniformly since this expansion can be done on the compact set A ri , 

and every geodesic variation is included in Mq. By the index lemma, we 

(jO)\ (t^) ~t 

have J-y n n (Jf i ) < J-y n n (£„ +1 ). Hence the desired inequality follows when 

o n = o. In the case o n ^ o, we have 

d (o) ,(o,^ a (tSi)) < d ( „) (o,o n ) + d (o) (o^X^i)), 



r/ 



(o,X a (t^)) = d ,, a) A0,0 n ) + d n( , a y(0 ni X a (t^)). 



Note that (tn , o, o n ) is uniformly away from Cutgx because of our choice 
of ro and Lemma 2.6. Therefore the conclusion follows by combining them 
with (3.8). □ 



Before turning into the next step, we show the following two comple- 
mentary lemmas (Lemmas 3.8 and 3.9) which provide a nice control of the 
second-order term A n in Lemma 3.7. Set A n = E[A n |^ n ,„i]. 

Lemma 3.8. Let (a n ) n6 N be a uniformly bounded ^ n -predictable pro- 
cess. Then 



lim a sup< 

a-)-0 



in probability. 



N+l 



j=n 



n, N e N, n < N < A/ (a) A 



o-R 



Proof. Note that the map (t,x,y) h-> Gt )X ,y(d(x,y)) is continuous on 
A ri . Since we have Gt )X ,y(d(x,y)) > on A ri , there exists K > such that 
K~ l < Gt :X ,y{d(x,y)) < K. This fact together with Corollary 3.6 yields |Aj| 
and |Aj| are uniformly bounded if j < ur. Since Y?j=i a jiAj ~~ A7) ls an =^n~ 
local martingale and <tr is J^-stopping time, the Doob inequality yields 



(3.9) lim a 2 sup 

Q ~ >0 0<N<N<,<*)/\a R 



N+l 



in probability. 
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U U' 

N=l n=0 



a 



N+l 
j=n+l 



>8 



c u u- 

N=l n=l 



a 



a sup 

0<N<NMA9 R 



I>j( a j-A 

JV+1 

^ajCAj-Aj) 









•i] 











jV+1 

^^■(Aj-Aj) 



> 



5 

>2 



Thus the conclusion follows from (3.9). □ 

Lemma 3.9. There exists a deterministic constant Co > being inde- 
pendent of a and R such that the following holds: 



A n +i < Co + - 



rd ia) (o,X^ n a) )) 



b(u) du, 



where b is what appeared in Assumption 1. 

Proof. By using (m + 2)E[(£ n , e*) (Cm e j)] = we obtain 



i=2 



(m-l)G' n (d(o n ,X«(t { n a) ))) 
G n (d( 0n „X-(4 a) ))) 



Note that we have 



(^(4 Q) ),7n) s(t W ) (X a (4 Q ))) - (Z(4 Q )),7n) 9(4 ( Q)) (On) 
d (a) (o„,X«(4 a) )) 

9(n ' a s (z(^)),7„) 9(t ( Q)) (7n( s )) 



du 



9(n > <v£ ) z(4«)), 7n ) 5(e)) ( 7n ( u ))^. 



Recall that, for (t,x,y) ^ CutsT, we have 

1 f d g(t)( x >y) 

dtd g{t) {x,y) = - 



(^5W)(7g(«).7g(«))^; 
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cf. [19], Remark 6. By combining them with Assumption 1, 
A n +i = dtd, t (c,)Jo,o n ) 

1 f d t.frh (°».* a (4 a) )) 

+ 2 1 d t g(t^)(j n (u),Uu))du 

+ (z(A%%) {a) (x-(m + (m " 1)G ; (d(o - x "j£ )))) 

} 2G n (d(o n ,,X-(4 a) ))) 

i f d (a) ( ,x«(4 a) )) 

(3.10) < - / 9( n ' 



2 . 

+ a t d 9(t ( Q)) ( , 0n ) + (^(4 a) ) i7n/ 9 ( t ( j "))l°nj 

+ 2j Q 3n RiC fl(t 0»))(7n(M),7n(«))<i« 

| (m-l)GUrf(on,A: a (4 a) ))) 
2G n (d(o n ,X a (d a) ))) 
Here we used the fact &(u) > in the case o n 7^ o. Note that 

/ • / \ ' / w i (m— l)G'(r) 

^C g{t ( n c) ) (jn(u),7n(u))du+^ ' wV 

is nonincreasing as a function of r. Indeed, we can easily verify it by taking 
a differentiation. Set 

Cf.= sup sup (\Z(t)\ g ( t )(x) 
te[Ti,T 2 ] xeB (t) {o) ^ 

+ sup (d t g(t)(V,V) + \mc g 

V£T X M 
\V\ Bit) <l 

By virtue of Lemma 2.2, C\ < 00 holds. By applying a usual comparison 
argument to G' n (ro)/G n (ro), we obtain 

d (a) K,x«(4 a) )) 

Ric g (4«) ) (7n («) , 7n 0) ) du 

| (m-l)G;(rf( 0n ,I"(fl a) ))) 

G n (d(o n ,X«(t^))) 

<Ci(r + coth(Cir )). 

Hence the conclusion with Co = Ci(l + 3ro/4 + coth(Ciro)/2) follows 
from (3.10). □ 
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In the next step, we will introduce a comparison process to give a control 
of the radial process. Let us define a function 99 on (2ro,oo) by 

<p(r) := C + - / b(u) du, 
z Jo 

where Co is as in Lemma 3.9. Let us define another function ip on (2ro,oo) 
so that ij) is a locally Lipschitz nonincreasing function satisfying ip{r) := 
2(r — 2ro) _1 for r G (2ro, 2ro + 1] and ip( r ) : = for r > 2tq + 2. Let us define 
a comparison process p a (t) taking values in [0,oo) inductively by 



/) a (Ti) := d g{Tl) (o,x ) +3r , 
# _ / ( Q ) 

p a (t) :=p a (t£>) + -^(«A n+1 + a 2 ( ¥ ,(P a (# )) + ^(p a (tW)))), 

The term is inserted to avoid a difficulty coming from the absence 

of the estimate in Lemma 3.7 on a neighborhood of o. By virtue of this 
extra term, p a (t) > 2tq holds for all t £ [Ti,T2] if a is sufficiently small. 
Let a' R and a' R be given by a' R := aji(p a ) and a' R := a~ 2 ([a' R \ a — T±) + 1. The 
following is a modification of an argument in the proof of [12], Theorem 3.5.3, 
into our discrete setting. 

Lemma 3.10. For 5 > 0, there exist a family of events (E?) a with 
lim Q _ s .oIP ) [-E'5 ] = 1 and 0, constant K{$) > with lim^o K(8) = such that, 
on Ef, 

d g{t) ( ,X a (t))<p a (t) + K(5) 
for t G [Ti,ctr A a' R A TJ[ and sufficiently small a relative to 5 and R~ l . 

Proof. It suffices to show the assertion in the case t = t { n a) for some 
n G No- Indeed, once we have shown it, Corollary 3.6(i) yields 

d g{t) ( ,X a (t))<e Ka \d g{[tia) ( ,X a ([t\ a )) + h(a)) 

<p<* tia+ K(5) + (e™ 2 -l)R + e™ 2 h(a) 

< pf + K(S) +a + (e Ka2 - l)R + e Ka " h(a) 

for t G [Ti,aR A T2]. Here we used the facts (p > and il>>0. From this 
estimate, we can easily deduce the conclusion. 



ti a) )) and p a (t^) 

by d n and p n , respectively, in the rest of this proof. Let us define a sequence 



For simplicity of notation, we denote d ,, {q k (o,X a {t { n>)) and p a {t { n>) 

9\Pn ) 
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of J^-stopping times Si by Sq := and 

S 2l+1 : = inf{j > S 2l \X a (& ] ) G Sif >} ( )} A tfW 



S 2i := inf{j > S 2l -i\X a (tf) i B^l(o)} A . 

Since p n > 2ro, it suffices to show the assertion in the case S21 < n < Sy+i A 
A a^j for some Z € No- Now Lemmas 3.7 and 3.9 imply 

d j+i ~ Pj+i < d j ~ Pj + a2 (f(dj) ~ f(pj)) + a 2 (A j+ i - Aj+i) + o{a 2 ) 
for j S [S2i,S2i+i Aa' R A a' R ). Here we used the fact tp > 0. Let f a be a C 2 - 
function on R satisfying: 

(i) /a|(— 00,— a) = 0j /a I (0,00) (■&) = x > 

(ii) / a is convex; 

(iii) a 2 sup^gK /^(x) = o(l). 

For example, a function / Q satisfying these conditions is constructed by 
setting 

f( x ) = / ftexpf- ° 2 jl(_ 1| i ) (s)dsdt, 
where a, b is chosen to satisfy 

/oo 6XP (~ 1^ ) ^ = 

bj j exp - _° g2 ^ (s) ds dt = 1 

and f a (x) := By the Taylor expansion with condition (iii) of f a , 

we have 

fa(d j+1 - p j+ i) < f a (dj - pj) 

(3.11) + a 2 f' a (dj - Pj )(<p(dj) - tpipj) + (A, - Aj)) 

+ o(a 2 ). 

Let C > be the Lipschitz constant of <p on [0, R]. Note that we have 

(3.12) f' a (dj - P j)(<p(dj) - <p(pj)) < C(d, - Pj ) + 

since <p is nondecreasing. Now by using (3.11) and (3.12) combined with the 
fact ds 2l — ps 2i < — a for sufficiently small a, we obtain 

(d n ~ Pn)+ < fa(d n ~ Pn) 

n—l n—1 

(3.13) < Ca 2 - Pj ) + + a 2 Y, f'M - Pj )(A j+1 - A j+1 ) 

j=S2k j=S2k 
+ 0(1). 
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Here the first inequality follows from condition (ii) of f a , and n < a^ 2 (T 2 — T\) 
is used to derive the error term o(l). Let Ef be an event defined by 

r k 1 

j=k 



Ef := I a 2 sup 

fc<fc'<A r ( Q )A 



<5 



Note that a,j = f' a (dj-i — Pj-i) is J^-predictable and uniformly bounded 
by 1. Thus, by combining Lemma 3.8 with (3.13), we obtain 

n-l 

(d n - Pn ) + < CO? ( d J ~ Pj)+ + 26 

j=S 2 i 

on E$ for sufficiently small a. Thus, by virtue of a discrete Gronwall in- 
equality (see [31], e.g.), 

(d n - Pn)+ < 25(1 + (1 + Ca 2 ) n ) < 25(1 + e c ^~ T ^). 
This estimate implies the conclusion. □ 

Corollary 3.11. For every R' < R, 

limsupP[(T fi < T 2 ] < limsupP[<r^ < T 2 ]. 

a— >0 a— »0 
Now we turn to the proof of our destination in this section. 

Proof of Proposition 3.4. By Corollary 3.11, the proof of Propo- 
sition 3.4 is reduced to estimate P[5^j < T 2 ]. To obtain a useful bound of 
it, we would like to apply the invariance principle for p a . However, there is 
a technical difficulty coming from the unboundedness of the drift term of p a . 
To avoid it, we introduce an auxiliary process p a in the sequel. 

Let (p be a bounded, globally Lipschitz function on R such that <p(r) = 
ip(r) + ip(r) for r 6 [2ro + JR — 1 ,i2]. Let us define an R-valued process p a (t) 
inductively by 

p a (Ti) :=d g{Tl) (o,x ) +3r , 

p a (t) :=p a (t^) + —^( a \ n+1 +a 2 v(p a (tt ) ))), te[t\t^]. 

We also define two diffusion processes p°(t) and p°(r) as solutions to the 
following SDEs: 

' dp°(t) = dB(t) + (<p(fP(t)) + V(P°W)) dt, 
p°(Ti) = d g (r 1 )(o,x ) + 3r , 

dp°(t)=dB(t) + <p(p°(t))dt, 
p°(Ti) = d g{Tl) (o,x )+3r , 
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where (-B(i))ie[Ti,T 2 ] i s a standard one-dimensional Brownian motion with 
B(Ti) = 0. We claim that p a converges in law to p° as a — > 0. Indeed, we 
can easily show the tightness of (p a ) a >o by modifying an argument for the 
invariance principle for i.i.d. sequences since (p is bounded. Then the claim 
follows from the same argument as we used in the proof of Theorem 3.1 
under Proposition 3.2, which is based on the Poisson subordination and the 
uniqueness of the martingale problem. 
Let us define rfR-.^i—t [T±, T 2 ] U {oo} by 

m {w) := inf{t G [T u T 2 )\w(t) < 2r + R- 1 }. 

Then we have 

n&R < t 2 \ < n<j R ( P a ) a m {p a ) < t 2 \ = na R {p a ) a m(p a ) < t 2 \. 

Since {w\or{w) A i]r(w) < T 2 } is closed in ^i, the Portmanteau theorem 
implies 

limsupPMp") A m (p a ) < T 2 ] < F[a R (p°) A m (p°) < T 2 ] 

= ¥[a R (p )A VR ( P )<T 2 }. 

Since p° is a diffusion process on (2ro,oo) which cannot reach the boundary 
by Assumption 1, the conclusion follows. □ 

3.2. Tightness of geodesic random walks. Recall that we have metrized 
the path space ^ by using d 9 (Ti)- To deal with the tightness of (^ a ) Q e(o,i) 
in ff, we show the following lemma, which provides a tightness criterion 
compatible with the time-dependent metric d g ^ . 

Lemma 3.12. (^ a ) a e(o,i) *s tight if 
1 



lim - lim sup sup 1 

<5-K) a~f0 neN 



sup d 9W (X Q (4 Q )),X Q ( S ))> e , 

A a) <s<(t' n a) +5)AT 2 



(3.14) 

a R = oo 

holds for every e > and R > 1 . 

Proof. By following a standard argument (e.g., [5], Theorems 7.3 and 7.4), 
we can easily show that (^ a ) a e(o,i) is tight if, for every e > 0, 

sup d g{Tl) (X a (t),X a (s))>e\ =0. 

L t<s<(i+<5)AT 2 J 



lim - lim sup sup ] 

<5->0 a ^0 te[T u T 2 ] 



Thus, by virtue of Proposition 3.4, (X a ) ae ^ ^ is tight if 
lim^-limsup sup P sup d g rr 1 )(X a (t),X a (s)) > e,a R = oo 

8^0 a ^o te[Ti,T 2 ] H<s<(t+8)AT 2 
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for every e > and R > 1. Given i? > 1, take Mo an d K as in Lemmas 2.4 
and 2.2, respectively. Then, for e < 1 and s,t G [Ti,T2], 

{d g{s) (X a (s),X a ([t\ a )) <e,a R = ^} 

C {d 9(Tl) (X a ( S ), X a (i)) < 2e K ^- T ^e,a R = oo}, 

if a is sufficiently small. Thus we have 

{ sup d g(Tl) (X a (t),X a (s))>e,a R = ^} 

4<s<(t+<5)AT 2 } 

c{ sup d g(s) (X a ([t\ a ),X a ( S ))>- 2 —±,a R = oo\ 

L[tJa<S<(L*Jc<+2(5)AT2 1 ) 

for a 2 < 5, and hence the conclusion follows. □ 

Proof of Proposition 3.2. Take R > 1. By virtue of Lemma 3.12, 
it suffices to show (3.14). Take Mq C M compact and k as in Lemmas 2.4 
and 2.2, respectively. By taking smaller e > 0, we may assume that e < fo/2, 
where fo = fo(Mo) is as in Lemma 2.6. Take n G No with n < . Let us 
define a ^-stopping time ( £ by 

C £ := M{k G N |n < A; < N^\d git M } (X a (t^),X a (tP)) > e}. 

Then, for sufficiently small a, 

{ sup d g(s) (X a (t^),X a (s))> 2e,& R = oo\ 

t^<s<(t^+S)AT 2 

(3-15) 

C {a (Cs - n) < 5,a R = 00}. 

Set p k := X a (t[ a) ) for fe G N and f(t,x) := d g{t) (p n ,x). Note that f 2 is 
smooth on {/ < e}. Let us define X' k by 

We claim that there exists a constant C > such that 

(3.16) /(4+Wi) 2 < /(4 a) >P*) 2 + 2a/(4 a) ,p fc )A' fc+ i + C« 2 

for k < Q A A r ^ Q ) on {<tr = 00}. Indeed, in the same way as we did to ob- 
tain (3.8), 

/(tfc+l.Pfc+i) 

< M°W + 2a/(4 a) ,p fe )A' fc+1 + a 2 (A^ +1 ) 2 

(3.17) 
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+ 2aV(4 a ^^)(^/(C.^) + (^(C)'7^ ; ) g( ^) ) (p fe )) 



+ a*f{tt\p k )r^ )) {J L+i ) + o{a< 



(4 a) ) 



Here o(a 2 ) is controlled uniformly. Let K\ > be a constant satisfying that 
the g(t)-sectional curvature on Mq is bounded below by —K\ for every t £ 
[Ti^Tz]. Such a constant exists since Mq is compact. Then a comparison 
argument implies 



TpnPfc 

Here the right-hand side is bounded uniformly if k < £ e A A^ a ' . The remain- 
ing estimate of the second-order term in (3.17) to show (3.16) is easy since 
we are on the event {<tr = oo}. Applying (3.16) repeatedly from k = n to 
k = Ce ) we obtain 



£ 2 <2a^/(4 Q) ,p fc )A' fc+ i + C<5 



k=n 



on {a 2 (C e - n) < <5,<5\r = oo}. Set A^ (a) := sup{fc G N |fc < aT 2 b + n}. By 
taking <5 < (2C) _1 £ 2 , we obtain 

{a 2 (C e -n) < 5, <7r = oo} 
(3.18) C fe/(4 a) ,p fc )A' fc+1 > f^,a 2 (C £ -n) <5,<T fl = ool 



. fc=n 



A? 



c I sup 5^/(4 a) »Pfc)i/ /(t ( 



.n<N<N^ a) k 



L {/(4 Q) ,P fc )< £ } A/fc+1> 4a^ r 



Set 



1 



We can easily see that \Y/.\ < 1 and SfcL n +i^fc ^ s ^A/-martingale. By [11], 
Theorem 1.6, with (3.18), we obtain 

F[a 2 (( £ -n)< S,a R = oo] 



< 



N+l 



sup 



> 



n<7V<A^ a) fc=n+l 



4oVm + 2 
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< 



cxp 



e 4 



8 v / m + 2(ae 2 + 4a 2 y/m + 2(Nf' - n)) 

-4 

< exp 



8 



8^/m + 2(ae 2 + 4y/m + 25) , 
Hence (3.14) follows by combining this estimate with (3.15). □ 

4. Coupling by reflection. For k £ K, let U a k be a one-dimensional Orn- 
stein-Uhlenbeck process defined as a solution to the following SDE: 

dU a , k (t) = -^U a>k (t)dt + 2dB(t), 

U a ,k{T\) = a. 

More explicitly, U a>k (t) = e -Ht-Ti)/2 a + 2 e k(s-t)/2 dj B( s ). H ere 5(t) is 
the standard one-dimensional Brownian motion as in the proof of Proposi- 
tion 3.4. 

Theorem 4.1. Suppose 

(4.1) 2(VZ(t)f + d t g(t) < Ric 9(t) +kg(t) 

holds for some fcGl. Then, for each x\,X2 E M, there exists a coupling 
X(i) := (X\(t), X2(t)) of two Jiff diffusion processes starting at (x\,X2) sat- 
isfying 

[t^T^^ > °] - ¥ [ T ^ <T U d g{Tl) (^)Mt) > 

rfg(Ti)(^l,^2) 



= A 

/or eac/i T € [Ti , T2] , where 



2^P{T - T 1 



V 27r ^-« [ t , fc = 0. 

In addition, for i = 1, 2, is a solution to the martingale problem associ- 

ated with the time-inhomogeneous generator Jzf^ and the filtration generated 
byX. 

Remark 4.2. (i) Our assumption (4.1) extends existing curvature as- 
sumptions in two respects. On the one hand, (4.1) is nothing but (1.1) when 
Z(t) = and k = 0. On the other hand, (4.1) can be regarded as a natural 
extension of a lower Ricci curvature bound by k. Indeed, Bakry-Emery's 
curvature-dimension condition CD(/c,oo) (see [2], e.g.), which is a natural 
extension of a lower Ricci curvature bound by k, appears in (4.1) when 
both Z(t) and g(t) are independent of t. 
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(ii) Given k > 0, a simple example satisfying (4.1) can be constructed by 
a scaling. Indeed, for a complete metric g whose Ricci curvature is nonneg- 
ative, g(t) = e~ k( - t ~ Tl ^g satisfies (4.1) when Z(t) = 0. 

(iii) Prom the first item in this remark, when Z(t) = 0, one may expect 
that (4.1) works as an analog of Bakry-Emery's CD(£;,iV) condition, which 
is equivalent to Ric 9 > k and dim M < N when g(t) is independent of t, 
instead of CD(k, oo) since dim M = m < oo in our case. However, the follow- 
ing observation suggests us that we should be careful: let us consider (4.1) 
in the case k > and Z{t) = 0. When dtg(t) = 0, the Bonnet-Myers theo- 
rem tells us that the diameter of M is bounded and hence M is compact. 
Moreover, the Bonnet-Myers theorem still holds under CD(k, N) in the time- 
homogeneous case; see [3, 4, 23]. However, when g(t) depends on t, it is no 
longer true that (4.1) implies the compactness of M. In fact, we can easily 
obtain a noncompact M enjoying (4.1) with k > for some g(t) by following 
the observation in the second item of this remark. 



By a standard argument, Theorem 4.1 implies the following estimate for 
a gradient of the diffusion semigroup: 

Corollary 4.3. Let ((X(t)) te [ TuT2 ], (F x ) xeM ) be a Jt? t - diffusion process 
with ¥ x [X(Ti) = x] = 1. For any bounded measurable function f on M , let 
us define Ptf by P t f(x) := K x [f (X (t))] . Then, under the same assumption 
as in Theorem 4-1, we have 



lim sup 

y->x 



P t f(x) - P t f(y) 



^2irp(t - Ti) z ,z'eM 



In particular, Ptf is dg^y globally Lipschitz continuous when f is bounded. 

Proof. Let X = (Xi,X 2 ) be a coupling of ^-diffusions (X(t),F x ) and 
(X(t),¥ y ) given in Theorem 4.1. Let r* be the coupling time of X, that is, 
r* :=inf{< G [7i,T 2 ]|X(i) G D(M)}. Let us define X* = (A7,X|) of (X(t),F x ) 
and (X(t),V y ) by 



X*(t) := 



X(i), ifr*>t, 
(Xi(t),Xi(t)), otherwise. 



Since r* is a stopping time with respect to the filtration generated by X, 
and Xi (i = 1,2) is a solution to the martingale problem associated with 
the same filtration, X* is again a coupling of ^-diffusion processes. Since 
{r* >T} = {inf Tl < t < T d g(t )(X(i)) > 0}, Theorem 4.1 yields 

p t f(x) - p t f(y) = nnxm) - f(x* 2 (t))} 

=E[(f(xi(t))-f(x* 2 (t))n^>t}] 
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<P[r*>t] sup \f( z )-f( z ')\ 
z,z'eM 

. / d g ( Tl ){x,y) \ 

\2y/P(t-Tl)J z,z'£M 

Hence the assertion holds by dividing the both sides of the above inequality 
by d g ( Tl ^(x,y) and by letting y—tx after that. □ 

As we did in the last section, let (^xy)x,yeM be a measurable family of 
unit-speed minimal <?(t)-geodesics such that jxy joins x and y. Without loss 
of generality, we may assume that jxy is symmetric, that is, 7^ (d g u\ (x, y) — 
s) = lyx{s) holds. Let us define rh X y :T y M —>T y M by 

m®v :=v- 2(v,j$) g (t)i$(dg( t )(x,y))- 
This is a reflection with respect to a hyperplane which is g(t)-perpendicular 
to 'jxy ■ Let us define mi] ■ T xM -»• T y M by 

mg(«):=mg((// (t ^)(d g(t) (x,y))). 

Clearly miy is a y(t)-isometry. As in the last section, let <J>W : M ->■ G^(M) 
be a measurable section of the y(i)-orthonormal frame bundle G^ (M) of M. 
Let us define two measurable maps $f 5 : M x M 0® (M) for i = 1, 2 by 



mg^P^.y), (:r,y)eMxM\L>(M), 
$W(x), (x, y) £ D(M). 



Take xi,X2 G M. By using $ 4 , we define a coupled geodesic random walk 
X Q (t) = (Af (t),X$(t)) by Af (Ti) = Xi and, for t E [t^ , , 

: = v / ^+2^ 4a)) (x«(4 Q )))e„+i, 

Af (i) := exp^ )} + « 2 ^(4 Q) ))) 

for i = 1,2. We can easily verify that Xf has the same law as X a with 

Xq =Xi. 

In what follows, we assume (4.1). We can easily verify that it implies 
Assumption 1. Thus, by Theorem 3.1, (X Q ) Q , > o is tight under Assumption 1. 
In addition, a subsequential limit X afe — > X = (Ai, X2) in law exists, and it is 
a coupling of two .^-diffusion processes starting at x\ and X2, respectively. 
We fix such a subsequence (ak)keN- In the rest of this paper, we use the 
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same symbol X a for the subsequence X° fc and the term "a — > 0" always 
means the subsequential limit "o^ — > 0." 

We will prove that the coupling X obtained as above is a desired one in 
Theorem 4.1. We first remark that we can easily verify that Xi (i = 1,2) 
is a solution to the martingale problem associated with the filtration gen- 
erated by X in the same way as in the proof of Theorem 3.1. Set a l R := 
an(dg(.)(o, Xf (•))) for i = 1,2. We fix R > 1 sufficiently large until the be- 
ginning of the proof of Theorem 4.1. Let Mq cMbea relatively compact 
open set satisfying (2.1) for 2R instead of R. We next show a difference 
inequality of d g ^(X. a (t)). To describe it, we will introduce several notation 

as in the last section. For simplicity, let us denote 7 n ,„, by 7-. 

p J _ 'xf(tl Q) )x 2 "(ti Q) ) J ,n 

Let us define a vector field V n+ \ along 7 n by 

K+l := //t" a)) (fn+l " (^l,7n> fl( 4»>)7„(0)). 
Take uel m Let us define A* +1 and A* +1 by 

r ._ / 2 (£l+lAn) g{t My if (2/1,2/2) £ D(M), 



A* 



2\/m + 2(^ n+ i,v), otherwise, 
d {Q) (x«(tl a) )) 



+ 4n' 1 \ V n+l) J 1 {X«(4 Q) )^D(A/)}- 

For 5 > 0, let us define t s : ^1 -»• pi,T 2 ] U {00} by 

7s(t«) := inf{t > Ti|u;(t) < 5}. 
We also define t$ by f<5 := r,5(d 3 (.)(X a (-))). 

Lemma 4.4. For n e No uwi/i n < iV^, we have 

^^)C*°(*&) < (1 + |)^ Q)/2 ^)(^(4 a) )) 

(4.2) +e^ Q) /2 (aA ; +1 + a 2 A:+i) 

+ o(a 2 ), 

w/iere n < fj A irj; A cj 2 ^ and a is sufficiently small. Moreover, we can control 
the error term o(a 2 ) uniformly in the position of~K a . 

PROOF. When (t&*%X a (£&^)) ^ CutsT, (4.2) is just a consequence of 
the second variational formula for the distance function combined with the 
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index lemma for Z^* n \ To include the case (ti a) ,X»(d Q) )) GCutsT and to 
obtain a uniform control of o(a 2 ), we extend this argument. Let us define H 
and pi,p 2 : H -> [T x , T 2 ] x M x M by 

# := {(t,x,y,z)\t G [Ti,T 2 ],x,y,z G M ,d g{t) {x,y) > 6, 

d g {t){x,y) = 2d g{ t ) (x,z) = 2d gi t ) {y,z)}, 

Pi{t,x,y,z) := 

p 2 {t,x,y,z) := (t,y,z). 

If q = (t,x,y,z) G if, then Pi(q),P2(q) ^ Cutsx since z is on a midpoint of 
a minimal g(i)-geodesic joining x and y. Since H is compact, pi(H) and 
P2(H) are also compact. Hence there is a constant r\ > such that 

inf {|t - t'| + d g(t) (x, x') + d g(t) (y, y') \ (t, x, y) G Pl (fT) U p 2 {H) , 

(t',x',y')eCutsT}>V- 
Take a > sufficiently small relative to rj and 5. Set 



Pn := In 



Pn '■= ex P 



2 



By the triangle inequality, we have 

d 9(t ( Q)) (x Q (4-))) = d ff(4a)) (xf(^)), Pn ) + d^^xfctW)), 
)(x a (t&)) < + d^dk*? (i a i)). 

Since (ti a \ X^i^),^) G H, we can apply the second variation formula 
to each term on the right-hand side of the above inequality. Hence we ob- 
tain (4.2). For a uniform control of the error term, we remark that 7 n is 

included in Mo, and the g(in^)-length of 7„ is bigger than 5. These facts 
follows from n < f g A A a 2 R and the choice of Mq . Thus the every calcu- 
lation of the second variation formula above is done on a compact subset of 
[Ti, T 2 ] x Mq x Mq which is uniformly away from CutsT- It yields the desired 
result. □ 

Let us define a continuous stochastic process U" on K starting at a by 
U«(t) := e~ kt / 2 a + ae"^ 2 ^ e^^A* + tZ±fLj& a) j 

As a final preparation of the proof of Theorem 4.1, we show the following 
comparison theorem for the distance process of coupled geodesic random 
walks. 
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Lemma 4.5. For each e > 0, there exists a family of events {E") a such 
that P[-£7"] converges to 1 as a — > and 

(4.3) <W X °W) < U2 g(Ti)0 , a{Tl)) (t) + e 

for all t £ [Ti , T2 A A a l R A cr^] on E° for sufficiently small a. 

Proof. In a similar way as in the proof of Lemma 3.10, we can complete 
the proof once we have found E" on which (4.3) holds when t = tn £ 

[Ti , T2 A i<5 A cr R A a R ] . Set A; +1 := E[A* n+1 \& n ]. Then £™ =1 e*^'- 1 /2 (A* - A* ) 

is an j^-local martingale. Indeed, A* +1 is bounded if tn < a R A a R and so 
is A* +1 . Let us define E° by 

( N+l 



*?H bu P £e^ /2 (A*-A*)<-^ 

jv<jvC«) 2a 2 (T 2 



In a similar way as in Lemma 3.8 or [16], Lemma 6, lim Q ,_ ! .oP[-E£ l ] = 1 
holds. On E£, we can replace a 2 e fc *™ a) / 2 A* +1 in (4.2) with a 2 e fc *"°°/ 2 A* +1 + 
e/(2(T2 — Ti)). Since we have (m + 2)E[(£j, e^)^, ej)] = 6ki, we obtain 

A^i<-^ (tfe0) (x«(^)). 



Thus an iteration of Lemma 4.4 implies (4.3) on E® when t = tn . □ 

Proof of Theorem 4.1. Take e G (0, 1) arbitrarily. Let R > 1 be suf- 
ficiently large so that 

limsupP[o-)j A <x| < T 2 ] < e. 

It is possible by Proposition 3.4. Set a := d g i Tt ){x\,X2). Take Te [Ti,T2], 
and let 5 > be 5 > 2e. Then Lemma 4.5 yields 

F[t s >T}< F[{t s >T}D E° n {a R A & 2 r > T}] + 2e 
<P[t 5/2 (C/-)>T] + 2 £ . 

Thus we obtain 

limsupP[f 5 >T}< limsupPf inf U?(t) > 5/2 

a^O a^O Ue[Ti,T] 

by letting e J, 0. Note that J7" converges in law to C/ a as a — > 0. Since 
{w G C([Ti,T 2 ] -> M x M)|r 5 (d ff( .)(w(-))) > T} 
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is open, and {w \ inf^m-jy w(t) > 5/2} is closed in C([0, T] 
manteau theorem yields 

inf d g(t) (X(t)) > s] < liminfPfo > T] 

Ti<i<T 9WV K " \ - a->0 L J 



the Port- 



< lim sup I 

a->0 



inf m(t)>5/2 
lte[T u T] a 



< 



inf UJt)>8/2 



Therefore the conclusion follows by letting 5 4- 0. □ 

5. Coupling by parallel transport. As a final part of the paper, we will 
see that we can also construct a coupling by parallel transport by follow- 
ing our manner. In the construction of the coupling by reflection, we used 



a map rrixy- By following the same argument after omitting rhxy in the def- 
inition of rrixy, we obtain a coupling by parallel transport. The difference of 
it from the coupling by reflection is the absence of the term corresponding 
to A* , which comes from the first variation of arc length. As a result, we can 
show the following; cf. [16]: 

Theorem 5.1. Assume (4-1)- For x\,X2 6 M, there is a coupling X(t) = 
(Xi(t), X2(t)) of two J£t- diffusion processes starting at xi and X2 at timeT\, 
respectively, such that 

d q(t) (Mt))<e- k{t - sy2 d q( AMs)) 



(t) 



"git) 

for T\ < s < t < T2 almost surely. 



It recovers a part of results studied in [1]. In particular, a contraction type 
estimate for Wasserstein distances under the heat flow follows. 



Proof OF Theorem 5.1. Let us construct a coupling by parallel trans- 
port of geodesic random walks X Q = (X^X^) starting at (x\,X2) 6 M x M 
by following the procedure stated just before Theorem 5.1. By taking a sub- 
sequence, we may assume that X" converges in law as a — > 0. We denote 
the limit by X = (X\,X2). In what follows, we prove 

Pf sup (e fc */ 2 d 5(t) (X(t)) - e fcs /X (s) (X(s))) > e] = 

'-Ti<s<t<T2 J 

for any e > 0. By virtue of the Portmanteau theorem together with Propo- 
sition 3.4, it suffices to show 



lim I 



(5.1) 



sup (e fc */X (t) (X Q (t)) -e ks / 2 d g(s) (X a (s))) > e, 

T!<s<t<T 2 



CO 
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for any R>1. We write d n := e ktn )//2 d . < a ) . (X a (i£^)) in this proof for sim- 
plicity of notation. For 5 > 0, let us define a sequence of ^ n -stopping times Si 
by So := and 

S 2 i+x -Mij^S^djKSjAN^, 

S21 := inf{j > S 2 i-i\dj > 25} A . 

Note that ds 2l _ 1 < 35 holds on {a R A a R = 00} for sufficiently small a. As 
mentioned just before Theorem 5.1, the difference inequality (4.2) holds with 
A* = when S^-i < n < S21 A A a R for some / € No- In this case, the error 
term o(a 2 ) is controlled uniformly also in I. Let us define an event Ef by 

Ef:={ sup £ e^ Q) /2(A*-A*)< A! 

Then, as in Lemmas 3.8 and 4.5, we can show ]xm ot - > o'P[E? ] = 1. On Ef n 
{ct}j A j| = 00}, we have d^ < d n + 5 for S2/-1 <n< N < S21 if a is suffi- 
ciently small. Moreover, for n < S21-1 < N < S21, 

dN ~ d n < (d N - d5 2i _i) + d s 2 i-i < 5( 5- 
In the case S21 < N < S21+1, we obtain d^ — d n < 25. Thus d^ — d n < 55 
holds for all n < A r on £^ n {<t]j A <r| = 00}. Take 5 > less than e/10. Then 
our observations yield (5.1) since d g ( t }(X. a (t)) — d g n t \\(X.(\t\ a )) becomes 
uniformly small on {a R A a R = 00} as a — > 0. □ 

REFERENCES 

[1] Arnaudon, M., Coulibaly, K. A. and Thalmaier, A. (2010). Horizontal diffusion 
in C7 1 -path space. In Seminaire de Probabilites XLIII. Lecture Notes in Math. 
2006 73-94. Springer, Berlin. 

[2] Bakry, D. (1997). On Sobolev and logarithmic Sobolev inequalities for Markov 
semigroups. In New Trends in Stochastic Analysis (Charingworth, 1994) 43-75. 
World Scientific, River Edge, NJ. MR1654503 

[3] Bakry, D. and Ledoux, M. (1996). Sobolev inequalities and Myers's diameter the- 
orem for an abstract Markov generator. Duke Math. J. 85 253-270. MR1412446 

[4] Bakry, D. and Qian, Z. (2005). Volume comparison theorems without Jacobi fields. 

In Current Trends in Potential Theory. Theta Ser. Adv. Math. 4 115-122. Theta, 
Bucharest. MR2243959 

[5] Billingsley, P. (1999). Convergence of Probability Measures, 2nd ed. Wiley, New 
York. MR1700749 

[6] Blum, G. (1984). A note on the central limit theorem for geodesic random walks. 

Bull. Aust. Math. Soc. 30 169-173. MR0759783 
[7] Chavel, I. (1993). Riemannian Geometry — a Modern Introduction. Cambridge 

Tracts in Mathematics 108. Cambridge Univ. Press, Cambridge. MR1271141 
[8] Coulibaly- Pasquier, K. A. (2011). Brownian motion with respect to time-changing 

Riemannian metrics, applications to Ricci flow. Ann. Inst. Henri Poincare 

Probab. Stat. To appear. Available at arXiv:0901.1999. 



TIME-INHOMOGENEOUS GEODESIC RANDOM WALKS 



35 



[9] Cranston, M. (1991). Gradient estimates on manifolds using coupling. J. Fund. 

Anal. 99 110-124. MR1120916 
[10] Ethier, S. N. and Kurtz, T. G. (1986). Markov Processes: Characterization and 

Convergence. Wiley, New York. MR0838085 
[11] Freedman, D. A. (1975). On tail probabilities for martingales. Ann. Probab. 3 100- 

118. MR0380971 

[12] Hsu, E. P. (2002). Stochastic Analysis on Manifolds. Graduate Studies in Mathe- 
matics 38. Amer. Math. Soc., Providence, RI. MR1882015 

[13] Ikeda, N. and Watanabe, S. (1989). Stochastic Differential Equations and Diffu- 
sion Processes, 2nd ed. North-Holland Mathematical Library 24. North-Holland, 
Amsterdam. MR1011252 

[14] J0RGENSEN, E. (1975). The central limit problem for geodesic random walks. 
Z. Wahrsch. Verw. Gebiete 32 1-64. MR0400422 

[15] Kendall, W. S. (1998). From stochastic parallel transport to harmonic maps. In 
New Directions in Dirichlet Forms. AMS/IP Studies in Advanced Mathematics 
8 49-115. Amer. Math. Soc, Providence, RI. MR1652279 

[16] Kuwada, K. (2010). Couplings of the Brownian motion via discrete approximation 
under lower Ricci curvature bounds. In Probabilistic Approach to Geometry. 
Advanced Studies in Pure Mathematics 57 273-292. Math. Soc. Japan, Tokyo. 
MR2648265 

[17] Kuwada, K. and Philipowski, R. (2011). Non-explosion of diffusion processes 
on manifolds with time-dependent metric. Math. Z. To appear. Available at 
arXiv:0910.1730. 

[18] Kuwada, K. and Philipowski, R. (2011). Coupling of Brownian motion and Perel- 
man's ^-functional. J. Fund. Anal. 260 2742-2766. 

[19] McCann, R. J. and Topping, P. M. (2010). Ricci flow, entropy and optimal trans- 
portation. Amer. J. Math. 132 711-730. MR2666905 

[20] Oshima, Y. (2004). Time-dependent Dirichlet forms and related stochastic calculus. 
Infin. Dimens. Anal. Quantum Probab. Relat. Top. 7 281-316. MR2066136 

[21] Philipowski, R. (2009). Coupling of diffusions on manifolds with time-dependent 
metric. Seminar talk at Universtitat Bonn. 

[22] Pinsky, M. A. (1976). Isotropic transport process on a Riemannian manifold. Trans. 
Amer. Math. Soc. 218 353-360. MR0402957 

[23] QlAN, Z. (1997). Estimates for weighted volumes and applications. Quart. J. Math. 
Oxford Ser. (2) 48 235-242. MR1458581 

[24] Stannat, W. (1999). The theory of generalized Dirichlet forms and its applications 
in analysis and stochastics. Mem. Amer. Math. Soc. 142 viii+101. MR1632609 

[25] Stroock, D. W. and Varadhan, S. R. S. (1979). Multidimensional Diffusion Pro- 
cesses. Grundlehren der Mathematischen Wissenschaften [Fundamental Princi- 
ples of Mathematical Sciences] 233. Springer, Berlin. MR0532498 

[26] Topping, P. (2009). £-optimal transportation for Ricci flow. J. Reine Angew. Math. 
636 93-122. MR2572247 

[27] VON Renesse, M.-K. (2004). Intrinsic coupling on Riemannian manifolds and poly- 
hedra. Electron. J. Probab. 9 411-435 (electronic). MR2080605 

[28] Wang, F. Y. (1994). Successful couplings of nondegenerate diffusion processes on 
compact manifolds. Acta Math. Sinica 37 116-121. MR1272513 

[29] Wang, F.-Y. (1997). On estimation of the logarithmic Sobolev constant and gra- 
dient estimates of heat semigroups. Probab. Theory Related Fields 108 87-101. 
MR1452551 



36 



K. KUWADA 



[30] Wang, F. Y. (2005). Functional Inequalities, Markov Semigroups, and Spectral The- 
ory. Mathematics Monograph Series 4. Science Press, Beijing, China. 

[31] Willett, D. and Wong, J. S. W. (1965). On the discrete analogues of some gen- 
eralizations of Gronwall's inequality. Monatsh. Math. 69 362-367. MR0185175 

[32] Zhang, Q. S. (2011). Sobolev Inequalities, Heat Kernels Under Ricci Flow, and the 
Poincare Conjecture. CRC Press, Boca Raton, FL. MR2676347 

Graduate School of Humanities and Sciences 

Ochanomizu University 

ohtsuka 2-1-1, bunkyo-ku 

Tokyo 112-8610 

Japan 

E-MAIL: kuwada.kazumasa@ocha.ac.jp 



